Converting genetic network oscillations into somite spatial pattern 
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In most vertebrate species, the body axis is generated by the formation of repeated transient 
structures called somites. This spatial periodicity in somitogenesis has been related to the tempo- 
rally sustained oscillations in certain mRNAs and their associated gene products in the cells forming 
the presomatic mesoderm. The mechanism underlying these oscillations have been identified as 
due to the delays involved in the synthesis of mRNA and translation into protein molecules [J. 
Lewis, Current Biol. 13, 1398 (2003)]. In addition, in the zebrafish embryo intercellular Notch 
signalling couples these oscillators and a longitudinal positional information signal in the form 
of an Fgf8 gradient exists that could be used to transform these coupled temporal oscillations 
into the observed spatial periodicity of somites. Here we consider a simple model based on this 
known biology and study its consequences for somitogenesis. Comparison is made with the known 
properties of somite formation in the zebrafish embryo . We also study the effects of localized Fgf8 
perturbations on somite patterning. 



PACS numbers: 05.40.-a, 05.60.-k, 73.63.-b 

I. INTRODUCTION 

Somites are transient structures that form a periodic 
growing pattern starting from the head (anterior) and ex- 
tending to the tail (posterior) of a developing vertebrate 
embryo which ultimately give rise to both the segmented 
vertebral column and to the musculature. Formation of 
somites is a rhythmic process characteristic of the species 
at given temperature. For instance, in the chicken em- 
bryo one pair of somites is formed every 90 min at 37 
°C and in the zebrafish one pair is formed roughly ev- 
ery 30 min at 28 °C. The origin of these intracellular 
oscillations have been identified as due to time delays 
in the transcription and translation of herl and her7 
genes in the case of zebrafish and the hes7 gene in the 
case of mouse. The total number of somites produced 
is conserved within a given species typically somewhere 
between 50 to 70 pairs of somites form on each side of the 
anterior-posterior (AP) axis, but can vary dramatically 
between species, thus the zebrafish develops 30 pairs of 
somites. As new somites bud from the anterior end of 
the presomatic mesoderm (PSM) that extends back to 
the primitive streak and tail bud, new cells are added 
at the posterior end by cell division from the tail bud 
keeping the size of the PSM constant as new somites are 
generated . 

Somitogenesis is one of the best studied process of pat- 
tern formation in the developing embryo, and a num- 
ber of models have been proposed to address the mech- 
anisms underlying the generation of such periodic pat- 
terns @, H MM S, H: the Clock and Wavefront model 
proposed by Coke and Zeeman in 1976 is perhaps best 
known, but others include Meinhardt's Reaction Diffu- 
sion model; and Stern's Cell Cycle model. All of them 
suppose that oscillations of genes and gene products oc- 



cur in the cells of the PSM from which the somites derive 
(Fig. [T]). The Clock and Wavefront model postulates the 
existence of a longitudinal positional information gradi- 
ent down the axis of vertebrate embryos. This gradi- 
ent interacts with the cellular oscillator stopping the os- 
cillations and producing a rapid change in locomotory 
and adhesive behavior of cells when they form somites 
[HQj]]. The oscillation in the PSM is the somite clock, 
and the moving interfaces at the anterior end of the PSM 
where the positional information reaches a critical value 
is called the wavefront. It is the interaction beween the 
clock phase and the positional information that controls 
somitogenesis. Here we suppose the positional informa- 
tion is supplied by the concentration gradient of the gene 
product Fgf8 released at the tail bud of the embryo and 
creating a linear gradient of morphogen along the AP 
axis @, Biologically the signal created by the Fgf8 
concentration interacts with the clock phase (the intra- 
cellular concentration of a gene protein) leading to cell 
differentiation and somite formation. At low Fgf8 con- 
centration, cell arrangement becomes more compact and 
the epithclialization process underlying somite formation 
begins. 

Meinhardt's Reaction Diffusion model also incorpo- 
rates the idea of an oscillatory mechanism in the PSM 
cells producing an alteration between cell states as in the 
Clock and Wavefront model. One of the main aims of 
this model is to account for somite structure itself which 
at even early stages after formation shows heterogeneities 
between the anterior and posterior portion of an individ- 
ual somite. To account for this structure the reaction 
diffusion model postulates that cells can be in one of two 
possible states: anterior a and posterior p and switch 
from one state to another until they reach a stable state 
in the presence of a morphogen gradient forming a spatial 
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pattern such as aaappaaapp..., where a repeated segment 
aaapp constitutes a somite. 

Finally Stern's Cell Cycle model assumes that the os- 
cillations in the cells arranged along the AP axis are in 
synchrony with each other, but cells in the anterior PSM 
are further advanced in phase than cells in the posterior 
PSM. Segmentation is hypothesised to occur when cells 
reach a certain time window in their cell cycle. The most 
important point raised by this model is the question of 
how cells are synchronized biologically. It appears that 
the answer to this question lies in the deltaC gene that 
codes for the Notch ligand at the heart of the intercellular 
Notch signalling pathway which when disrupted destroys 
cell cycle synchrony All these models raise important 
questions, but recent experiments on how Fgf8 signalling 
controls somite boundary position and regulates the seg- 
mentation clock appears most in accord with the 
clock and wavefront model. 

Based in the Clock and Wavefront model combined 
with key ingredients of the known biology, we have there- 
fore constructed a deterministic one dimensional model 
and compare the resulting spatial patterning with the 
somites for zebrafish embryo. We use in our model the 
recent proposal for the mechanism producing the intra- 
cellular oscillations as due to delays in the synthesis of 
mRNA focussing specifically on the zebrafish genes herl 
and herl and their gene products the protein molecules 
Herl and Her7 

The paper is outlined as follows: Section [II] is devoted 
to the description of the proposed model and contains 
the results derived from it. The transcription in the 
synthesis of mRNA is essentially a noisy process and is 
considered for haploid cells in the section ITTT1 In section 
TVl we study the effect of a local perturbation of the 
Fgf8 wavefront on the spatial pattern of somites. Finally 
a conclusion and discussion is presented inlVl 



II. DETERMINISTIC ONE DIMENSIONAL 
MODEL 

The model consists initially of a linear array of N star t 
cells that represents the starting size of the PSM with 
intracellular oscillators driven by the herl and herl 
genes and their associated gene product proteins . The 
oscillations between nearest-neighbor cells are coupled 
and synchronized via Notch signalling [l^, The 
cell-cell communication produced by means of the 
deltaC gene product Notch crossing the cell membrane 
wall. 

In each cell in the array the gene mRNA concentrations 
m k for the genes k = herl, herl, deltaC and their asso- 
ciated gene products the protein concentrations pk obey 
the sets of coupled kinetic equations (see [l|): 
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FIG. 1: An schematic illustration of somite formation within 
the Clock and Wavefront model. Gene and gene product con- 
centrations in the cells of PSM experience temporal oscilla- 
tions. In the top part of the diagram, the positional informa- 
tion wavefront is illustrated together with the position of the 
determination front which we assume represents the position 
at which the Fgf8 concentration gradient falls below a critical 
value. The somites form in an anterior to posterior order as 
the wavefront advances from head to tail. 
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Here the symbol i denotes the cell position in the linear 
array; i' goes from 1 to nn the number of near neighbor 
cells of i (nn can be 1 if i is at an end or otherwise 2); T mk 
is the delay time from initiation of transcription to export 
of the mature mRNA m k into the cytosol; T Pk is the delay 
between the initiation of translation and the emergence 
of the complete protein molecule p k ; a k represents the 
protein p k synthesis rate per mRNA molecule; b k is the 
rate of protein p k degradation; c k is the rate of mRNA 
m k degradation; while the function f k represents the rate 
of production of new mRNA molecules m k that is given 
by the Michaelis-Menten kinetics 
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where we used the notation if = t — T mk ; and 
<j)k(i,t) = p k (i,t)/p c " t . The parameters rO k , rd k , rh k , 
and rhd k that add up to 1 , represent the relative weights 
of k transcription that is unregulated, regulated by 



deltaC protein alone, regulated by her protein alone, 
and regulated by her and deltaC proteins, respectively. 
p^ %t is the critical number of molecules of Herl or 
Her7 protein per cell, for inhibition of transcription 
if k is herl or herl respectively, while it represents 
the critical number of Notch molecules required for 
activation if k is deltaC. The mRNA concentrations for 
the kind different of genes decrease when the k protein 
concentration are larger than their critical values p c k Tlt 
if k is Herl or Her7 or when the deltaC protein Notch 
concentration is below its critical value p^itaC- ^hc 
concentrations of intercellular signalling molecules in 
the cells are coupled between the nearest-neighbour nn 
cells by means of the function /& through the deltaC 
protein Notch concentration of their neighboring cells 
Pdeitac{i', t - T m J), as it can be seen in the Eq. ©. 
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FIG. 2: Concentrations of mRNA molecules of gene herl 
and its associated gene product protein Herl as a func- 
tion of time for any cell of the linear array. Sustained os- 
cillations of both concentrations are obtained for the fol- 
lowing values of parameters involve in Eqs. JTj and 
a k = 4.5, b k = 0.23, c k = 0.23, K k = 33, r0 k = rd k = 0, 
with k = herl, herl, deltaC ', rhheri,her7 = 0, rh de itaC = 
1, rhd her i ther7 = 1, rhd her i >her7 = 1, rhd de i ta c = 0, 

Ph7rl,her7 = 4 0, PdMa = 1000, T mhcrl = 12, Tm hcr7 = 7.1, 

Tm deltaC = 16, T Phcrl = 2.8, T Phcr7 = 1.7, T VdMaC = 20 min. 
(Inset) The oscillations are stopped when the corresponding 
cell is reached by the somitogenesis wavefront. The concen- 
trations of protein and mRNA molecules of genes herl and 
deltaC have a similar behavior as the concentrations of rriheri 
and Pheri shown in this figure. 

The concentrations of both protein and mRNA 
molecules, pk and m k , are set equal to zero initially in all 
Nstart cells of the linear array. The temporal evolution 
of pk and rrifc are calculated with Eqs. (JlJ. Sustained 
oscillations of mRNA and proteins are obtained for 
parameters and time delays T„ lk and T Pk close to the 
zebrafish estimated ones as it is shown in Fig. [2] [l[. 
The mRNA and protein oscillations can occur without 
interactions between neighbor cells, because the delays 
that take part in feedback loops are intracellular. The 
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FIG. 3: Concentration of mRNA molecules of gene herl as 
a function of cell number in the somites. The concentrations 
were calculated using Eqs. (fl| and ([2]) and using the same 
parameter values as in the Fig. [2] (see caption). The con- 
centration of rriheri belonging to cells that were reached by 
the wavefront were frozen. In our simulation the wavefront 
advances at a speed of one cell every two minutes. We see 
the simulation at time t — 900 minutes at the end of somito- 
genesis in the zebrafish. At t — 900 minutes 29 somites have 
formed. All the other concentrations of protein and mRNA 
molecules in the cells (namely p k with k=Herl, Her7, Notch 
and rrik with k = herl, deltaC) have the same qualitative 
behaviour as the concentration of rriheri- 



interactions through f k connecting neighbor cells are, 
however, crucial as they result in synchronous in phase 
oscillations in the linear array representing the PSM as 
a whole to occur. 

We suppose at time t = all the cells in the starting 
PSM to be oscillating in phase, and therefore we begin 
with the linear array of N star t cells and integrate them 
for the time necessary to obtain coherent synchronous 
oscillations. In this paper N star t = 280, and in its syn- 
chronous oscillating state it represents the time t = in 
our model. The growth of the PSM is simulated in our 
model by the adding of new cells to the initial group of 
cells at one end of the linear array. This growth veloc- 
ity Vg row th is a variable parameter in our model, and we 
chose to add cells at a rate of one cell every 5 minutes 
in agreement with measured growth rates [lOj. Thus, it 
generates sequentially along the antero-posterior axis the 
necessary cells for forming the future segments of the em- 
bryo: the somites that are created at the anterior portion 
of the PSM, see Fig. [TJ In our simulations once the initial 
PSM has been obtained as a group of cells oscillating in 
phase, the caudal motion of the wavefront starts. Again 
the wavefront velocity v waV e is a parameter in the model 
that controls the rate at which cells reached by the wave- 
front stop oscillating and thus freezing the gene and gene 
product concentrations at their values when they meet 
the wavefront (see inset of Fig. [5]). The observed spa- 
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tial shift in the phases of different cells is thus obtained 
by the wavefront catching different phases when passing 
through different cells, a biological implementation of the 
clock and wavefront model. 

What value should we choose for v wave ? The simplest 
argument would suggest v wave = v growth- This supposes 
that the FGF creates a time indepedent distribution on 
a fast timescale compared to growth, and therefore the 
wavefront is simply reading off a critical value of the 
FGF concentration at which the intracellular oscillator 
is frozen. If we chose this value, however, we would not 
create around 30 somites in 900 minutes as we would 
expect [10| . Indeed, it is important to note that the 
regularity of the somite pattern and the size of somites 
is very sensitive to changes in N start and v wave - For 
example, if N star t = 20 and v wave = v growth k, 5 min 
per cell, the number of somites is about 10 for 900 
minutes. Here we report instead on our simulations with 
v-wave — 2 minutes per cell. Then when the simulation 
finishes in the time t — 900 minutes we find that 29 
somite pairs have been created. A wavefront advancing 
at a rate of 1 cell/2 minutes creates a spatial wavelength 
A = T ce uv wave « 15 cells in length. What is the process 
of cell differentiation after the cell meets the wavefront 
now depends on the biology. Here to be specific we have 
assumed that when the rriheri level is above its average 
value a somite is formed. Averages of the concentrations 
of proteins and mRNA of the frozen cells are calculated 
at this final time. With this choice a somite is roughly 
6 — 7 cells in length. On the other hand the neighbor 
cells that have protein and mRNA concentrations lower 
than average form the space between somites. With this 
choice the remaining 8 — 9 cells form an intersomite gap. 
Of course, it is also possible to interpret the simulation 
as forming somites of length 15 cells with a differentiated 
structure between the anterior and posterior portions of 
each somite. Now as the simulation proceeds the PSM 
reduces in size. Indeed, for the parameters chosen here 
N s tart — 280 cells, but N en d = lOcells, because while 
180 cells are added in the 900 minutes of the simulation 
at the posterior end of the PSM, a total of 450 cells 
differentiate into new somites, It is not surprising 
therefore that phase coherence is lost towards the end of 
somitogenesis, and the spatial array of somites obtained, 
though approximately periodic, decrease in size towards 
the tail (Fig. EJ). 



III. STOCHASTIC MODEL 

Our simulations above assumed that no noise was 
present in the cells during transcription and translation, 
but as Lewis has pointed out, the transcription step in- 
volved in protein synthesis is essentially a stochastic pro- 
cess because of the small number of molecules involved 
in the cell [T^: A DNA molecule can randomly have a 
Herl, Her7 or Delta Notch dimer bound to its regula- 



tory site. When such a dimer does bind, transcription 
is forbidden; whereas if no protein is bound at a regu- 
lation site of the DNA molecule, transcription occurs at 
the possible maximal rate. 

The stochastic model for haploid cells is then con- 
structed by taking the deterministic one dimensional 
model and considering the possible states (bound and 
dissociated) of all regulatory sites of the herl, herl and 
deltaC genes as stochastic variables [lj. A DNA regu- 
latory site can have a Herl, Her7 or Delta Notch dimer 
bound or not (i. e. either transcription forbidden or 
maximally free of repression). A gene regulatory site 
then makes stochastic transitions between the bound and 
unbound states with a certain probability in any time in- 
terval. 

The protein and mRNA concentrations are given by 
the Eqs. {1} as in the deterministic case, but the function 
fk now describe a stochastic process and are therefore 
modified to: 

fk(pheri(i,t'),p her7 (i,t'),p de ita(i',t')) = K k {r0 k 

+rd k [l-Zl D) (i',t')}+rh k tl H) k(i,t') 

+rhd kV i H \i,t')[l-r 1 i D) {i\t')}} 

(3) 

where the parameters r0 k , rd k , rh kl and rhd k were de- 
fined previously, below of Eq. @. The functions £, k D \ 

£*^> Wft*^ an d r n k D \ with k = herl, herl, deltaC, are 
random variables (rv) taking a value 1 when no protein 
is bound at the regulatory site and when it is bound, 
with probabilities that depend on protein concentrations 
in the cell i and in its neighbor cell i', at time t' — t — T mk 
when synthesis of m k (t) begins. These probabilities at a 
given time t are thus conditioned by the values the ran- 
dom variables take at an earlier time t' . Specifically the 
probability that any random variable (rv) takes value 1 
at the time t + At depends on the state of the appropriate 
regulatory site at time t (bound or unbound). 

p[%(t + At) = l-P rv (t + At)\(P rv (t)=0), if regulatory 

site is bound at time t 
plKt + At) = P rv (t + At)\(P rv (t) = 1), if regulatory 

site is unbound at time t 



In Eq. [¥] p™p represents the probability that any ran- 
dom variable rv has value a at t+At if it had value /3 at t 
(rv can be <^ ff) , rf^ or n^, and a = l,/3 = 0, 1); 
while P rv (t) is the probability that the regulatory site 
is free of repression. This probability obeys the master 
equation 

dP rv /dt = k r v ff (l-P rv )-k r v n P rv . (5) 

These kinetic coefficients are discussed in [l[ (Supplemen- 
tal Data) , where fc™ the rate constant for protein binding 
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is a MIchaelis-Menten function of the protein concentra- 
tions, while fc^Jy is the rate constant for dissociation of 
the repressor proteins from their regulatory DNA bind- 
ing sites (typically biologically fc^J * rj 1 min^ 1 implying 
a mean lifetime of about 1 min for the repressor bound 
state). 

Solving the above differential equation results in 
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the forms P rv (t + At) -- 
P rv (t)exp[-v rv At], where v r 



t™(l - exp[-vAt]) 
(K v ff + Kn) and u™ 



k r V ff/v rv . Inserting these solutions in Eq. Q one obtains 



p\l{t + At) = u rv {\ - e [ - vr " At] ), if regulatory site is 

hound at time t 
plUt + At) = u™(l-e [ - t, "" Atl ) + e [ - t,r " Atl ! 

site is unbound at time t 



— 1 1 1 r 

t . 

m herl in each cell of lineal array 
somites 




if regulatory 



100 200 300 

number of cells 



(a) 



200 



U50 



e 

■B ioo 



e 

c 
U 



50 



r 

without noise 
low noise level 
high noise level 




200 

time [min] 



400 



FIG. 4: Concentration of messager RNA molecules of gene 
herl as a function of time for any cell of the linear ar- 
ray. Comparison of results obtained using the determinis- 
tic one dimensional model described in section [TT] (solid line) 
and the stochastic model for a low noise level (dashed line, 
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= 10 and k o d ^ taC = 1) and for a high noise 

V {H) v <D l AH) 

level (dotted line, kJ^T 117 — kj}}' 1 ^ = 1 an d k^£l taC = .1). 
The parameter values are the same that in Fig. [2] The con- 
centrations of protein and message molecules have the same 
behavior as the concentration of rriheri shown in this figure. 



In the same way, the probability that any random vari- 
able (rv) takes value at the time t + At is given by 



Pll{t + At) 
p r l{t + At) 



(1 - u rv ) (l - e 1- ""'*' 1 ) , if regulatory 

site is unbound at time t 
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FIG. 5: Concentration of message molecules of gene herl as 
a function of the number of cell at the linear array. The 
concentrations were calculated using the stochastic model for 

„( H ) „(£>) 

haploid cells with a noise level given by kJi e T x = & D // rl/7 — 
Ah) 

5 and fc o ^ ,taC = .5, and assuming the rest of the parameters 
equal to the used ones in the Fig. [5] The wavefront advances 
at a speed of 2 minutes by cell and 26 somites are formed at 
t = 900 min, in the simulation. All the other concentrations 
of protein and message molecules in the cells {i.e. pt with 
k — herl, herl ,deltaC and m*, with k = herl ,deltaC) have 
similar behavior as the concentration of ra^eri- 



Expressions ^ and joined with the Eqs. |1| and 
^ allow to calculate the concentration of protein and 
mRNA molecules at anytime instant for the linear array 
of cells considering transcription as a stochastic process. 
There is a correspondence between this system with noise 
and the deterministic system: statistical averages of the 

random variables Cfe £> \ £fc > Vk an< ^ ^fc^ m the limit 
of large k r Q v ^ and fc™ tend toward the behavior of the de- 
terministic system. In Fig. [51 the concentration of m^ri 
versus time calculated using the deterministic model (de- 
tailed in the section [TTJ) and the model with noise for two 
different levels of noise are compared. A random variabil- 
ity in the amplitude and shape of individual oscillation 
peaks can be seen for both results obtained with noise. 
In addition, as the amplitude of noise increases the oscil- 
lations of m^eri are not in phase with the deterministic 
model oscillations. 

Finally, a somite pattern for the stochastic model is 
shown in Fig. [5l The wavefront advances at a speed 
of 2 minutes per cell freezing the oscillations of protein 
and mRNA concentrations and the simulation finishes at 
time t=900 minutes. Thus, a pattern of 26 somites is 
obtained with the same constant values specified in the 
caption of Fig. [2] and also used in Fig. [3) It is possible 
to increase the number of somites decreasing the delays 
within the range of allowed values (see [l|). 



IV. PERTURBATION OF THE 
SOMITOGENESIS WAVEFRONT 

The concentration of Fgf8 signalling is high in the cau- 
dal PSM and drops between the intermediate and the 
rostral PSM as indicated by the top part of the Fig. [1] 
In this paper we assume that somite formation can begin 
when the Fgf8 signalling falls below a critical level, and 
this level moves towards the caudal end of the PSM with 
a constant velocity, the somitogenesis wavefront v wave . 
Note we assume this velocity is not the growth rate of 
the PSM but in fact is faster v wave > v growt h, we shall, 
however, assume the same velocity in both the determin- 
istic and stochastic models. Now it is known that a tran- 
sient manipulation of the wavefront in zebrafish embryos 
alters the size of the somites [TtJ HH : larger somites re- 
sult when there is transient inhibition of Fgf8 signalling, 
whereas smaller somites result in the presence of tran- 
sient activation. Chemical inhibitor and transplantation 
of Fgf8 beads are used for altering the level of Fgf8 sig- 
nalling that regulates the position of the wavefront within 
the PSM [H 111 . 
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Assuming that such perturbations in Fgf8 concentra- 
tion directly result in changes in the wavefront velocity, 
we will modify the wavefront velocity in our model and 
analyze their consequences. The perturbed pattern of the 
resulting somites depend on how the wavefront velocity 
Vwave changes. In Figure examples of such perturba- 
tions are shown that have been obtained with the deter- 
ministic one dimensional model by altering the wavefront 
velocity (Fig. [5]). In each case the other parameters in 
the model required for finding the protein and mRNA 
concentrations are the same as in Fig. [3] Specifically 
as in the unperturbed system, the wavefront initially ad- 
vances a rate of 2 min. per cell but between the cells 
100 and 160 the wavefront is perturbed. The velocity is 
increased to 1 min per cell in the cell interval (100, 160) 
and larger somites are formed (see Fig. |6ja)) in com- 
parison with the results of Fig. [3] When the velocity is 
decreased to 4 min per cell in the same interval of cells 
smaller somites are formed (Fig. Eljb)). In general al- 
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FIG. 6: Concentration of gene herl mRNA as a function of 
the number of cell in the linear array. The concentrations 
were calculated using the deterministic version of the one di- 
mensional model and assuming the same parameter values as 
in the Fig. [2] (see caption) and|3j The wavefront advances at 
velocity of 2 min per cell in both figures, but between the cells 
100 and 160 there is a perturbation in the Fgf8 levels resulting 
in a change in the wavefront velocity: the velocity decreases 
to 1 min per cell in (a) and increase to 4 min per cell in b). 
The patterns of somites are altered by the perturbations (see 
Fig. EJ. 



terations of the wavefront velocity perturb the resulting 
pattern somite in a grossly irregular manner. 



V. CONCLUSIONS 

We have carried out simulations (see Figure [5|) us- 
ing the deterministic version of our model that result 
in a growing approximately spatially periodic sequence 
of somites by combining the clock and wavefront model 
with the temporal oscillation ofgene transcription and 
translation suggested by Lewis [Ij for the zebrafish em- 
bryo (Fig. [3j. We also simulated the transient manip- 
ulation of the wavefront in zebrafish embryos (ToL ITU . 
Larger somites are formed when the wavefront velocity is 
locally increased whereas smaller somites results if there 
is a local decrease in the wavefront velocity (Fig. [6]) . In- 
deed, we believe that these are the first simulations that 
show in a biologically plausible manner the local varia- 
tion of somite size due to external perturbation of the 
Fgf8 concentration. 

We also noted in section IIII1 the consequences that 
gene regulation is in reality a noisy process which is 
likely to be crucial in the real developmental situation 
because of the small number of intracellular molecules 
involved. Our simulations also raise a number of 
questions for the biology os somitogenesis. First, as 
gene regulation is noisy it is important to study the 
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consequences of this noise for somite patterning. Second, 
the interplay between growth at the caudal end of the 
PSM and the velocity of the somite development at the 
rostral end is crucial for the number, size and shape 
of the resulting somites. It would be very useful there 
to have more biological information of the relationship 
between these processes. Finally, we assume that the 
Fgf8 concentration acts as a morphogen creating an 
intracellular signalling cascade that ultimately controls 
transcription and translation of the gene products. 



We do not have detailed information of this process 
and simply assume that the Fgf8 acts as an on/off 
switch. It is certainly possible, however, that a more 
complex dynamics may be involved and it would be 
useful to have more information on Fgf8 signalling in 
this context. In conclusion, further developments in the 
model will depend on knowledge of this biology and on 
the pertinence to somite formation and differention on 
the phase of the oscillation cycle [H . 
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